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Abstract: We discuss the black hole effective action and define its static subsector. We 
determine the induced gravito-static polarization constants (electric Love numbers) of static 
black holes (Schwarzschild) in an arbitrary dimension, namely the induced mass multipole 
as a result of an external gravitational field. We demonstrate that in 4d these constants 
vanish thereby settling a disagreement in the literature. Yet in higher dimensions these 
constants are non- vanishing, thereby disproving (at least in d > 4) speculations that black 
holes have no effective couplings beyond the point particle action. In particular, when 
l/{d — 3) is half integral these constants demonstrate a (classical) renormalization fiow con- 
sistent with the divergences of the effective field theory. In some other cases the constants 
are negative indicating a novel non-spherical instability. The theory of hyper geometric 
functions plays a central role. 
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1 Introduction 

The word "stereotype" carries a negative connotation for judging an individual based on 
a group to which he or she belongs rather than basing on actual individual traits. ^ Yet, 
in this work we shall describe a useful treatment of compact objects including black holes 
which can be considered as stereotyping (see figure 1). 

Most of this paper is devoted to a determination of the coefficients of induced gravito- 
static polarization, known as the electric Love numbers [2], ^ of a black hole in 4d and in 
higher dimensions. Consider a non-trivial gravito-static space-time and consider placing a 
small non-spinning black hole at an equilibrium point. As a result of the tidal gravitational 
force the black hole will deform, mass multipole moments will be measurable from far 
away, and could be interpreted as changes in the internal mass distribution of the black 
hole. To leading order the induced multipole moment is proportional to the external field 
multipole at the equilibrium point and we shall refer to the constant of proportionality as 
the coefficient of induced gravito-static polarization, or Love number. 

^ Stereotype (dictionary definition [1]) - A preconceived and oversimplified idea of tlie cfiaracteristics 
which typify a person, situation, etc. 

■^Named after the British mathematician A. E. H. Love who defined them first in 1911 in the context of 
Newtonian gravity. 
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STEREOTYPES 

Exist for a reason 



Figure 1. 



As a concrete example, consider two very massive and approximately static bodies (see 
figure 2) . On the line joining them there is a point of gravito-static equilibrium. We place 
a small BH at that point and wish to measure the induced mass multipole moments from 
far away. Another concrete example, one which is exactly static, is a caged black hole, 
namely a small black hole within a larger extra dimension (see for example [3], [4], [5] and 
reference therein) . In this case the black hole reacts to a gravitational field created by itself 
(through its "mirror images"). 

These constants are basic linear response coefficients for black hole physics, much like 
compressibility in elasticity. Apart for the intrinsic interest, knowing these constants would 
have several applications. First, they are needed for solving the binary motion. While for 
4d black holes such effects appear first (in the two body effective action) only at order 5PN, 
they are much more relevant for other types of compact objects and moreover in very high 
dimension the leading effect for BHs appears already at order 1"'"PN, namely only a little 
over IPN. Another possible application is the study of extended black objects in higher 
dimensions and their instabilities, as we shall find out. 

Even though black hole Love numbers could have been defined following Schwarzschild's 
discovery of the black hole metric [6] it was only rather recently that works on the subject 
appeared. [7] computed the second Love number for various models of neutron stars and 
discussed the possibility of measuring it using gravitational wave detectors. Then in June 
2009 three works appeared [8, 9, 10] which addressed the black hole Love numbers, yet two 
issues remained open 

• While both [8, 9] independently found that a certain calculation yields zero contribu- 
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Figure 2. Experimental set-up for measuring the induced gravito-static polarization constants. 
The black hole (gray) is positioned at an equilibrium point between two fixed stars (black). Its 
mass multipole moments are measured on a far away envelope (dashed sphere). 



tion to the Love numbers, they disagreed on whether another contribution should be 
included, one arising from the effective field theory side and which was too difficult 
to compute at the time. ^ 

• [8] made an interesting speculation that a black hole may have no non-minimal world- 
line couplings whatsoever. ^ 

We shall consider the determination of the Love numbers from the point of view of 
the effective field theory approach to GR [11, 12] including the ideas of [4].^ We started 
working on this topic over two years ago but were held back by the above-mentioned 
apparent ambiguity in their definition. 

The strategy of this paper is to calculate the Love numbers for an arbitrary space-time 
dimension following the idea of "dimension as a parameter of GR" , see [14] and references 
therein including [15]. We shall discover an interesting dependence which will also shed 
new light on the 4d issues. 

We start in section 2 by discussing the concept of the black hole effective action and 
defining its static subsector. We define the induced polarization constants. In section 3 
we describe a concrete physical set-up for its measurement together with the correspond- 

^[8] wrote "The question of computing the 'correct' value of ki [the Love numbers] for a black hole is a 
technically much harder issue which involves investigating in detail the many divergent diagrams that enter 
the computation of interacting point masses at the 5-loop (or 5PN) level." while [9] "boldly proclaimed" 
that the tidal Love numbers of a black hole must be zero. Finally [10] wrote "there are subtleties inherent 
in any definition of the multipole moments of BHs, so that there is currently no unambiguous determination 
of the ki Love number of BHs." 

^"This vanishing suggests, but in no way proves, that the effective action describing the gravitational 
interactions of black holes may not need to be augmented by nonminimal worldline couplings." [8] 

^See [13] for early precursors of the EFT approach to GR. 
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Figure 3. The idea beliind tlie black liole effective action: the full black hole space-time geometry 
is replaced by a point particle together with some effective interactions with its slowly varying 
background. 

ing computation. In section 4 we perform the "microscopic" calculation, namely we find 
the appropriate non-asymptotically-flat deformation of the Schwarzschild solution and we 
determine the Love numbers. The results are interesting and we proceed to discuss them 
in subsection 4.2. In section 5 we discuss divergences and their physical meaning including 
counter-terms, classical RG flow and cancellation of divergences in the EFT. In section 6 
we look at the calculation from an EFT perspective and reproduce several results with 
Feynman diagrams. We end with a summary and discussion in section 7. 

2 The black hole effective action 

Consider a soliton of size ro moving in a background with typical length scale L ^ tq. Given 
the hierarchy of scales it is natural to use the point particle approximation. However, this 
approximation is certainly limited, as it does not retain almost any of the object's physical 
properties, namely the finite-size effects. To improve upon that one uses either a matched 
asymptotic expansion between a near zone and a far zone, or an effective field theory 
(EFT) approach which is the more modern tool (see for example [4] which followed the 
EFT approach [11]). The idea is illustrated in figure 3. 

A black hole (BH) is a soliton of interest in Einstein's theory of gravity. In [11] a 
matching procedure was outlined which produces the black hole effective action. Later [4] 
defined it essentially® as follows 

SfuU,eff[x'',e'^;Qint;g] = SEFT,eff[x'',e'^;Qint;g] (2.1) 

:= Sbuik {gBH{x'^,e'X;Qint;9)) 

■= Sbuikig] + SBH[x'',e'^;Qint;g] 

This definition means that the effective theory is built such that its effective action for a 
point particle in a given background, is identical with the same effective action in the full 

®[4] did not list the internal degrees of freedom as they appear here - those were later discussed in 
[16, 17]. 



- 4 - 



theory with a fuU black hole metric inserted into the same background. Accordingly, given a 
space-time metric g, one must solve for a black hole metric qbh which describes a black hole 
(namely a horizon) embedded on an appropriate world-line given by x^(r) in g and having 
the prescribed spin and internal degrees of freedom, e'^,Qint respectively. The full bulk 
action (Einstein-Hilbert -|- gauge fixing) is evaluated on gsH- The resulting expression has 
a (standard) bulk contribution and a contribution localized on the world-line. The black 
hole effective action is defined to be the latter. 

We now wish to define the static sector of Sbh- In this sector the background metric 
is static (namely time independent and time reversal symmetric), the black hole location 
is fixed and it cannot rotate nor excite its internal degrees of freedom. Therefore in this 
sector 

Sbulk (gSHig)) = SfuU,eff[g\ = SEFT,eff[g] = Sbulkig] + SBH,st[g] (2-2) 

and the BH effective action depends only on the background (but not on any collective 
or internal coordinates). For self-consistency the location of the black hole must be an 
equilibrium point of the gravitational background, namely a static geodesic, which without 
loss of generality we take to be the origin, 0. 

In the static limit it is convenient to use the Non-Relativistic Gravitational (NRG) 
field decomposition [4, 18], namely a temporal KK reduction, which replaces the Einstein 
field gfj_u by three fields: the Newtonian potential the gravito-magnetic vector potential 
Ai and the spatial metric jij according to 

ds^ = e^^idt - Ai dx' f - e-'^'t'l^ jij dx'dx^ (2.3) 
where we conveniently denote 



d:=d-3. (2.4) 



In the static limit Ai = 0. 

We define a Lagrangian by 



SbhM'^^^] = j dt CBH,st[^,l] ■ (2.5) 

We expand J~-BH,st around flat space-time ^ij = 5ij,4) = 0. In general the leading term in 
Sbh is the universal point-particle action 

Sbh = -m J dr + . . . 

dr := ^ g^i,{x)dxi' dx" (2.6) 

Therefore in the static sector 

^%,st = (2.7) 

Perturbations around fiat space-time are parameterized by (j), aij where 7jj = 5ij +aij. 
Terms proportional to the equations of motion can be removed by a field re-definition. In 
particular 

dii^, a,(^|o (2.8) 
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are redundant terms, the first being a linearized equation of motion in the bulk (assuming 
a vacuum space-time) and the second being the linearized geodesic equation (no gravita- 
tional field at 0). We now wish to examine the leading terms beyond the point particle 
approximation, namely the finite-size corrections. There can be no terms linear in a due 
to the spherical symmetry (of the Schwarzschild black hole) - for instance a term linear in 
(j) would represent a mass multipole for Schwarzschild. 

The most general terms quadratic in (f) can be written as follows 



K ^ 



Cbhm = -m + -Y,h\d^'<l^\ +••• (2.9) 

1=2 

where Ii := (ii, . . . ,i/) is a multi-index and each i runs over the spatial directions i = 
1, . . . , ci — 1 and K will be defined shortly below in (2.12). We start with (j) (rather than 
(Tjj) because the leading finite-size effect for the post-Newtonian two-body effective action 
comes from the quadrupole (/ = 2) term above. ^ The sums starts with 1 = 2: / = is 
absent since 4> is only derivatively coupled, and Z = 1 is absent due to (2.8). We refer to the 
constants A/ as the induced gravito-static polarization constants, or tidal Love numbers. 
Their normalization here was chosen in a rather natural way, such that (for generic d, I) 

Qi^ = -Xi<j)j^ <Pj^ := di,cP\o (2.11) 

where Q/, is the mass multipole at the origin induced by the background field multipole 

In particular, the factor K is chosen such that the kinetic term for (p is S D —K/2 J di(f)di(f). 

From(4.8) we have 

The sign convention is such that in weak (Newtonian) gravity A > 0. Most of this paper is 
devoted to computing A/ for all d. The dimensions of A/ are 

= M L^^ = L^^'^^ . (2.13) 

Returning now to our title "black hole stereotyping" we see that according to effective 
field theory it is useful to "stereotype" a BH when it is viewed from afar by a series of 
properties each with a specific order. ^ The mass is most relevant, followed by the spin 
which is non-static and cannot be seen in (2.9). The next property and the first finite-size 
one is A2, and so on ... 



^In addition in 4d Uij cannot fluctuate since in 3d space Ricci flatness implies a flat space. 

* Despite considerable effort so far we were not able to find a clear enough definition of the standard 
normalization of the tidal Love numbers to allow comparison. 

^Actually there is some freedom in how we count the order. Apart from derivative counting we may also 
assign dimensions to the fields t/i, o". For instance in PN should be considered order m while a is order 



-6- 



3 Measuring Love 



Here we translate the general definition of the Love numbers A (2.9) into an experimental 
measurement set-up and a corresponding computation. 

In the experimental set-up, following the description in the introduction, one positions 
a non-rotating black hole at an equilibrium point of a weak gravito-static field (such that 
the typical distance for non- linear effects L^i is much larger than tq). The Newtonian 
potential (p is taken to be of the form of a pure /-multipole, namely (j) \s a, homogeneous 
function of degree I in the spatial coordinates and satisfies the Laplace equation ducj) = 0. 
One then measures the induced mass /-multipole from the asymptotic form of cj), and 
computes A as the constant of proportionality in (2.11). 

The computation starts with a background field (p{x) satisfying the flat space eq. 
of motion dii<j) = 0. As usual we may perform a decomposition into spherical harmonics 
(j) = Ylim 4'lmif)Ylm{^) where in a general dimension m stands for a multi- index. A general 
solution to the radial equation in flat space is of the form (pi = Ar^ + Br~^~'^. Imposing 
origin regularity we have 

(pl = ^rK (3.1) 

Now we should obtain (pBH^ the field in the presence of a black hole which asymptotes 
to (p and is regular on the horizon (for clarity we suppress now the index I). For r <^ L^i 
we can write (pBH as a linear combination 

(pBHir) = ^[Mr) + U2ir)) (3.2) 

where 



rM 1 + . . . O 



-d 



l + ...0(^)) (3.3) 



and A is some constant. 

Considering the equivalent definition of A (2.11), and recognizing that the induced 
mass multipole at the origin Qi is proportional to (j)X, the coefficient of <j)2, while (pj is 
proportional to (j) we find that 

A = iVA (3.4) 
where the normalization constant is found to be (see appendix A) 

N = An — ^ . (3.5) 

2'r(f + /) 

We would like to make two comments about our procedure (3.2-3.5). 
Possible ambiguity. In the flat space theory (pi, (p2 in (3.3) have no corrections. However, 
in the presence of corrections the definition of (p\ becomes ambiguous to adding a multiple 



Again at distances r such that ro r <C L^i 
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of 02 whenever the powers in the two series (3.3) overlap, namely whenever (2/ + d)/d G 
N = 1, 2, 3, . . . . We may express this as 




where we defined 



Though in general / need not be integral, in 4d / is nothing but /, and this ambiguity is 
essentially the one noted in [10]. 

A simple and effective solution is to observe that for generic values of d no such issue 
exists, and then the special values can be treated as a limit. Alternatively, one must use 
the general definition (2.2) to separate the contributions to A originating from two different 
EFT vertices (2.9): the m vertex and the A vertex. 

We note that in EFT language the condition (3.6) means nothing but that A has a 
dimension which is an integral power of m. 

Later we shall find that the small parameter for the asymptotic expansion (3.3) is 
actually m^/r^'^, at least in the background of a free scalar. Hence the condition for 
possible ambiguity (3.6) becomes more restrictive, namely 



Consistency with earlier terms of the EFT. The general definition of the static 
effective BH action (2.2) requires that the EFT action reproduces the effective action of 
the full theory. In the present context it means that corrections to (at least those which 
dominate (^2) should be reproduced in the EFT by the m- vertex alone. Indeed the equation 
of motion involves the black hole metric which can be fully reproduced from the m- vertex 
while the A- vertex reproduces the boundary condition (horizon regularity). In section 6 we 
shall see several examples for this consistency. 

4 Microscopic computation 

In our first computation we consider a minimally coupled scalar field rather than the 
Newtonian potential (j). We do that for several reasons (i) this case is simpler (ii) the 
gravitational case will turn out to be similar and (iii) this case actually describes the 
master field for tensor mode perturbations [19], namely perturbations of the spatial metric 
in d > 4. 

The action is 



(3.8) 




We note that = 1 instead of (2.12). We wish to solve for ip in the background oi a, d 
dimensional black hole in Schwarzschild coordinates 




(4.2) 
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where tq is the Schwarzschild radius. 

It is convenient to change coordinates r — )• X where 



(4.3) 



namely r^X = 1 — f }^ 

Decomposing into (normahzed) spherical harmonics the action becomes^^ 



S 



d a 



dX 



(4.4) 



where / was defined in (3.7). The equation of motion is 







dxfdx 



X^ 



(4.5) 



We note that the equation of motion depends on d,l only through the combination /. 
Moreover, this static equation is analytic in X while the time dependent equation contains 
the term X"^/*^ = which is non-analytic for d > 5. It would be interesting to explain 
these simplifications. 

Solution. The equation of motion (4.5) has three singularities all of which are regular: at 
the singularity, the horizon and asymptotically. Therefore the solutions can be expressed 
by the hypergeometric function. From hereon we use units such that tq = 1. The solution 
which is regular at the horizon is 



iPiX) = X-'^ Fi-l,-l,l; 1- X) = F{1 + 1,-1,1; -{1- X)/X) = Pf 



X 



X 



(4.6) 



The second equality uses the Pfaff identity (B.4) and are the Legendre polynomials. 

In order to read A we must expand around r = oo, namely X = 0. We do that using 
(B.5) and we use the Gamma function identities (B.7) to simplify the expression, finally 
arriving at 

1 r2 (/"+!) 



A 



tan7r/r2'+'^ 



(4.7) 



where we restored units using powers of tq. 



Isotropic coordinates. At first sight dimensional analysis (3.6) allows for divergences 
when 2/ + 1 = 1,2,..., yet these occur only for 2Z + 1 = 2, 4, . . . . This suggests that the 
small parameter can be taken to be rather than m. 



^^We record the standard relation between the Schwarzschild radius ro and the mass m, namely rg = 
16irGm/ 

^^Note that we consider the static limit, therefore the integral over time decouples and we suppress it in 
what follows. 
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Indeed, this is the case. The static Einstein-Hilbert action, when expressed in terms 
of NRG fields (2.3), reads [4] 



s[4>,iij\ 



1 



and enjoys the following parity symmetry 



1 + ) f^ix)^^<^)^j(P + R[-f] 



7*j ^ 7ij ■ 



(4.8) 



(4.9) 



When supplemented by m — >• — jtt, it is a symmetry also of the relevant interaction term 
—mcf). In gauges which respect this symmetry, such as the isotropic and the harmonic 
gauges we may conclude that = 7jj(m^) while (p{—m) = —(p{m). Moreover, Tp couples 
only to 'jij as can be seen by expressing its action (4.1) in the static case in terms of NRG 
fields 

~d'^xY^{x)dii;dj^p (4.10) 



Therefore is indeed the small parameter for ijj, thereby justifying the improved dimen- 
sional condition (3.8). 

The Schwarzschild gauge which we used here (4.2) provides a simple metric but does 
not respect this symmetry and therefore this parity property was not manifest. It turns out 
that there is a different representation of the black hole metric which preserves the parity 
symmetry while providing equations of motion which are quite similar to Schwarzschild 
coordinates. These are the isotropic coordinates 



ds' = {fjufdt'-ff{dp' + p'dn') 

/± := 1 ± (po/p)' 



(4.11) 



This metric is written in NRG form and satisfies the parity property (4.9). 

The coordinate transformation which takes us from isotropic to Schwarzschild is given 

by 



1 + 



(4.12) 



In particular, the location of the horizon in isotropic coordinates is at po which is related 
to ro through 



The static action in isotropic coordinates now becomes 

d 



S 



where 



z 



\ozi^i\ +— Fzl 



(4.13) 



(4.14) 



(4.15) 



-10- 



The solution to the equation of motion which is regular at the horizon is given by^^ 

ip{Z) = Z"'>(l/2, -1, 1; 1 - Z^) . (4.16) 
In terms of pQ our result (4.7) simplifies a bit and becomes 

r2(f+i) 



A,; 



r(f+i r / + ! 



, f d{2l+l) 

tanvrt 



(4.17) 



Combining with (3.4,3.5) this completes the computation of the induced polarization 
constants for the case of a scalar field background. 

4.1 Gravitational polarization 

Our starting point is the master equation for gravitational perturbations of scalar type in 
d dimensions [20] which are the d dimensional generalization of the Zerilli equation [22] 
and describe perturbations associated with the Newtonian potential. In our notation the 
equation reads 



d^ ^ 2{d-l)X + 2 d ^ {d-l)X -{l-l){l + d+l) 



Y{X) = 



(4.18) 



dX^ ' dX{X-l) dX ■ d'^X^{l-X) 
where we continue to suppress the index I, namely Y = Yi. We change variables into 

Y:=rY (4.19) 

and the equation reads 



fdjc - 2dx 



X2 



Y = . 



(4.20) 



This equation can be gotten from the following action 



S[Y] 



dX 



dxY 



X2 



Y 



(4.21) 



As in the case of a minimally coupled scalar the equation is analytic in X and depends 
on l,d only through /. We comment that this action encodes very economically the rather 
longer Zerilli potential, and that its form suggests that Y is the time component of a 
contravariant vector, namely Y = for some vector V. 
The solution is 



Y{X) = X~^ F{l-I,-l,2; 1 -X) = F(/ + l,-f,2; -{1-X)/X) 



Expanding around X = we find 



r(or(/ + 2) ~ rf(2f+i) 

tanvr/ Pq 



rff+i r [+1 



1 + y ) 



(4.22) 



(4.23) 
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Induced Gravitostatic polarization constants 

A 




Figure 4. A graph of the raw gravito-static Love numbers A of a black hole as a function of 
I = l/{d— 3) where I is the spherical harmonic index, and d is the total space-time dimension. Units 
are such that po ~ 1 where po is the location of the horizon in isotropic coordinates (4.13) 



where po was defined in (4.13). 

The graph of A is given in figure 4. 

In order to relate the result for A (4.23) with our original definition (2.9) we must find 
the relation between Y and the Newtonian potential (p. When going to the EFT we must 
take the flat space limit of Y. Since there is only one gauge invariant scalar the limit of 
Y and 4> must be proportional to each other, namely Y — ?■ g{r)(j) for some g{r). Next 
we observe that in the flat space limit the equation of motion for Y becomes = Al". 
Comparing with the (flat space) equation of motion of (/>, namely = Ac/) we conclude 
that g{r) = const and actually we may normalize our definition of Y such that g{r) = 1. 
Altogether we found that in flat space 

Y ^(j) (4.24) 

(it would be nice to show this directly from the definition of 1", see also [23, 24]). Therefore 
we can proceed and combine (4.23) with the normalization factor (3.4,3.5) to obtain our 
main result - the induced gravito-static polarization constants or a Schwarzschild black 
hole. 

4.2 Discussion 

We can now resolve the issues concerning 4d Love numbers which were listed in the intro- 
duction: 

• 4d Love numbers (of electric type) vanish. 

^^This transformation is known in the mathematics literature as a quadratic transformation of the hy- 
pergeometric function, see for example [21]. 



- 12 - 



Indeed in 4d d = 1, / is always integral and tanvr/ = 0. This settles a disagreement 
in the literature (described in the introduction), validating the intuition of [9]. We 
note that in general [8] were also right to consider an additional EFT contribution, 
and indeed we shall find it to be non-zero for half integral I. Yet in 4d it happens to 
vanish. 

Our approach was able to avoid the danger of ambiguity in definition since for generic 
dimensions none exists, and the 4d limit is smooth. Thus it can be considered a 
dimensional regularization. The 4d limit also serves as a first test of our results. 

The vanishing of Love numbers is quite surprising especially that the black hole 
horizon is known to deform, see for example [8, 10, 25, 26]. It implies a certain 
infinite rigidity of the mass distribution of a 4d black at least with respect to linearized 
gravito-static perturbations. 

• Black holes do have non-minimal world-line couplings for d > 4 (where I could be 
non- integral) . This settles negatively the speculation made in the abstract of [8], at 
least for d > 4. 

The tangent factor in (4.23) has several implications 

• A vanishes whenever I is integral. 

• A has a pole whenever I is half-integral. This is related to a classical renormalization 
fiow and will be discussed below. A second test of our result will come from a 
comparison of the residue of such poles against EFT divergences. 

• A(/) oscillates and changes sign. In particular it can be negative. The interpretation 
of a negative value is not clear at the moment. A negative compressibility of an 
elastic material implies an instability. However, a black hole cannot be divided and a 
negative specific heat, for example, does not imply an instability of a black hole, but 
rather of the black string [27]. This instability is spherisymmetric but we believe the 
during pinching this symmetry would be spontaneously broken. We speculate that 
negative Love numbers are related to this last instability. 

• In 5d there are only zeros and poles. 
We make some additional observations: 

• The factor in (4.23) which includes the Gamma functions has quite a smooth behavior 
for I > 0, for large I it approaches 1 and corrections in l/Z can be computed 



r2(/ + i) 




3 3 




(4.25) 



r(/ + i)r([+| 



) 



d 32P 



This behavior suggests that a theory can be formulated in the / — )■ oo limit. 
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• At integral values of I the solutions are especially simple ~ they are polynomials in 
X, since when either o or 6 are 0, —1, —2, . . . the hypergeometric series F{a, b, c, ; x) 
is finite. A useful explicit formula is 

F(c + k, b, c; X) = ^.li^X^-^X'^+'^-i (1 - X)-' (4.26) 

for A; = 0,1,2,... 



5 Diverging Love 

In this section we examine more closely the case where I is half integral and both Love 
numbers (4.17,4.23) diverge. The situation is familiar in quantum field theory where it 
leads to RG flow, and here it provides another example, quite pleasing, for classical RG 
flow see [11, 28, 29]. 

Cancellation of divergences. While the gravitational A diverges (4.23), Y, the horizon 
regular wave-function, (4.22) remains finite, of course. Looking at the asymptotic expansion 
according to (B.5) and substituting a = 1 — l,b = —I, c = 2 for Y (4.22) we have 

F(l - [, 2; 1 - X) = r(2)r(2/ + l) _ ^ _^ 

r(/ + i)r(/ + 2) 

+ lHliz^[zilx2'+iF(f+l,f+2,2f+2;X) (5.1) 

r(i-or(-0 

The divergence in A originates from a divergence in the factor r(— 2Z — 1) in the second 
summand. However, at the same time the c parameter in the first hypergeometric function 
becomes a negative integer ci = —21 — ?• —n which implies a divergence in that function. 
Quite generally the diverging part in such a case is given by 

F{a,b,-n + €;X) = (5.2) 

- "°r(:)r(6ir('!;t'2r + ^) + " + M + n + l. „ + 2; x, + o (.°) 

where n = 1,2, .. . and T{—n + e) = (— )"/(n!e) + O (e*^). Using this relation with a = 
1 — l,b = —I, n = 21 one confirms that the two divergences indeed cancel as e — )• 0. 

Log terms. As a byproduct of the cancellation of divergences and the presence of di- 
mensional constants a log(r) term is generated in Y{r). Indeed, as the divergences cancel 
at order O (e~^) the power law X'^^^^ = X""*"^^*^ in the second summand (reaction term) 
generates a finite O (e*^) logX term multiplied by the residue of the pole. Hence the coef- 
ficient of the log term, which is interpreted as the beta function, is directly related to the 
coefficient of divergence (the pole residue in our dimensional regularization). 

The same phenomenon is seen also from a solution of the differential equation of 
motion (4.20) through the substitution of an asymptotic power series Y = '^^Y^X^ . For 
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half integral I the resulting recursion relation at order k = 2/ + 1 becomes • ^2[+i ~ 
source ^ 0, which is known from the mathematical theory of differential equations to 
signal the appearance of log terms. 

Finite, RG-fiowing part. Next we compute the finite part. To do that, we must separate 
A into a counter-term and a finite part. The counter-term must cancel the EFT divergence. 
The finite part depends of course on a choice of renormalization scheme, and we use the 
standard MS (Minimal Subtraction) scheme. As we saw log terms get generated and 
accordingly the finite part depends also on the scale in which it is evaluated. So altogether 
we must keep both the pole residue and the finite part in A, but not higher terms in the e 
expansion. 

For convenience we turn to the case of a free scalar. Let us now assume that / = n/2 is 
half integral, that is n = 1,3, 5, .... Then the induced polarization (4.17) diverges. However, 
according to (4.16) the wave function tp is finite. As will be explained later, the divergent 
part of induced polarization is cancelled by another divergence inherent to the full solution 
(4.16). This means that A^ is given by the sum of an infinite part and a finite term which 
we now turn to compute. To separate these contributions, we introduce the following 
definitions 

m = L^nir , e = d — do , (5.3) 

where L is an arbitrary length scale and has length dimension - do. Of course, m must 
be L-independent, therefore rrir fiows with L according to 

dm drrir 
L— = L—^ = -emr. (5.4) 

UiIj cLLj 

Substituting these definitions into (4.17) and expanding in 6, yields"^^ 



A, 



4, r(H±i)r(n±3) 

2njo _ ^ n\ndo - 2) - n ^ _ 
e (do + l)(n + l) 



2HnL + do logvr — 2do log(L//3i 



'0, 



+ do n{n + 1) ^{do/2 + 1) - ^^(do(n + l)/2) + 0{e) ) , (5.5) 



where Ha is the a-th harmonic number and ^{z) is the digamma function. Notice that we 
should keep L'' in front of the above expression to maintain correct dimensions for A^ in 
{d + e)-dimensional spacetime. In the special case / = 2, do = 4 (n = 1), we obtain 



A„ 



l=2,do=4 



11^2 7 \ , L ^ 
- + - log - + — - 7 + log — + 0{e] 
e 2 \ TT bU J po 



(5.6) 



The pole in e corresponds to the divergence. In the EFT approach this term represents 
the so-called counter term which has to be introduced into the theory in order to make the 



^*Note that po depends on d, therefore one has to express it in terms of m first and only then proceed 
with expansion in e. 
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observables, e.g. the wave function ■0, finite. In contrast, the finite piece in the expression 
for represents the renormahzed value of A^, i.e. 



Afdoin+l) p2 (n 

A7(L)-^^^o 1 2 



47r r(^)r(^) 

— n „ f 

2Hn + do logvr - 2do log(L/po) 



n^{ndQ — 2) — n 



X ( 2n^Hn-i - 2 — -■ ■■ n 

+ do n{n + 1) ^-(^0/2 + 1) - '^^{do{n + l)/2)) . (5.7) 

In particular, while A^ is certainly L-independent, the induced polarization A™(-L) exhibits 
a classical RG flow 

- ^^L~ - r(^)r(^)^° ■ ^ 

6 EFT side 

In this section we demonstrate and confirm the equivalence of the EFT description by 
mirroring some of the previous matching calculations on the EFT side. 

For convenience we consider the case of a free scalar field (4.1) and we choose the 
isotropic coordinates (4.11). 

Integrating out short scales of order po, leaves us with the following effective action 



^ oo 

SEFT,eff{(r, (p, Tp) = Sbulk{<y-, (p, Tp) + SBH,st{(y-, 4') + l^^^k \d^'"p\ , (6.1) 

where^^ 



2 

1=2 



Sbuik{<y, ^) = S[^] + \ ' ' / dx^^'a— {day / dx'^^'a-^ (50)^ , (6.2) 

and the indices are contracted with flat metric biy To simplify the gravity action (4.8), we 
set 7jj = abij for the background metric in isotropic coordinates and used the following 
transformation rule for the [d — l)-dimensional curvature scalar 

= e^'^bij R[-f] = -e-2'^[2(d + 1)V^lo + d{d + l)diiod'uj] , (6.3) 

where oj = log(cr/2). Note that the full action is quadratic in ip and symmetric to constant 
shifts, ijj ^ ip + const. This explains why only quadratic terms in the derivatives of ip are 
present in the above effective action. 

In what follows we match the quadrupole A2. Introducing the following perturbation 

^^^ + ip i) = ipijx'x^ , Qi = (6-4) 



^For simplicity we set G = 1 for the rest of this section. 
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p 

/WWWV 



1 

>VWV/ 



52tt 1 



-k 



(d-2)(d-3) p' 



. P 

>VW\/ 



d-3 



-(d-3)(2TT)'-' lij'iK. 5(k+p) 



- m 



A,ilj''p^Pj 



Figure 5. All necessary diagrams to match A2 in isotropic coordinates. Solid, dashed and wavy 
lines represent propagators of 0, ip and a respectively. Cross indicates insertion of tp. Wave numbers 
flow into the vertex. 



Y = Y + 

^ EFT ^ ^ 



.+ 




(a) 



(b) 




+ . 



Figure 6. Diagrams contributing to the asymptotic value of the scalar field -0. 



into the effective action (6.1), we evaluate the asymptotic expansion of the scalar field, 
tpEFT, and subsequently match the result with the corresponding expansion of the full 
solution (4.16). The Feynman rules can be obtained by expanding (6.1) in the weak field 
approximation, i.e., a = o" — 1, ^ 1- All necessary rules for our needs are listed in 
figure 5. In what follows we use dimensional regularization to evaluate various diagrams. 
All necessary formulas to accomplish the computations can be found, e.g. in the appendices 
of [28]. 

The contribution of the induced quadrupole A2 to ipEFT is shown in figure 6(a) and is 
given by 



Fig.6(a) = 4^43^V'. 



(6.5) 



pd+4 

For d > 5 this contribution dominates all other possible Feynman graphs, since according 



-17- 



to (2.13) in this case A2 ~ m"' with 1 < q < 2, whereas as explained above the smah 
parameter of the theory is m^. Expanding (4.16) in Z ^ 1, yields 



. ^ ^ f 1 + V" '-^^^ ^ tan(2./d) + (6.6) 

V vpy r(3/2 + 2/d)r(i/2 + 2/(i) y 



Matching with the EFT result, we obtain 

r((i/2 + 2)r(3/2 + 2/(i)r(l/2 + 2/d) ^ ' ^'^0 
This result agrees with (4.17) upon substituting / = 2 and taking into account (3.4). 

The divergence at I = 2, d = 7. Let us now consider d = 4, namely Z = 2 in 7d.^^ In 
this case we should include the diagram in figure 6(b), since A2 ~ m? 

Fig^6(b) = v;f«)"nM^. (6.8) 
^ ' \pj r(3 - rf/2) ^ ' 

This diagram diverges in d = 4. Thus in the language of quantum field theory, we have 
to renormalize the indefinite (bare) parameter A2 such that becomes finite. We now 
demonstrate this procedure. 

Using (5.3) to expand figure 6(b) in e, yields 

Fig.6(b) = - 1 (^^y + ^ - 87 - 8 log vr - 16 log £ + 0(e)) . (6.9) 



where 7 is the Euler constant. To cancel the pole in e, we redefine A2 as follows 

e 



I A A 

A2 = L^pI ( Xf' + — + -|- + other possible poles in e I , (6.10) 



(i) 

where X2 , i > 1 are e-independent constants adjusted to eliminate all poles in e which 
appear in ip, while A2''^ depends on an arbitrary scale L and should be matched to agree 
with the full solution in this case. This is the so called minimal subtraction (MS) scheme. 

Pictorially, the above redefinition can be depicted as in figure 7. Namely, one can 
think about it as decomposition of the bare parameter A2 into a finite piece proportional 
to A2'''' and an infinite part associated with Aj^^ In the QFT literature the latter is called 
a counter term. 

Substituting (6.10) into (6.5) and expanding in e, yields 



Fig.6(a) = I (^) + 6Af + IIaW - 3A^^^^ - 3A^^^ logvr - 6A^^^ log ^ + 0{e) 



(6.11) 



"The more general case oi I — 1/2 with arbitrary d is quite similar. 
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(a) 



(b) 



Figure 7. Pictorial representation of decomposition of (6.10). (a) represents the finite part, whereas 
(b) corresponds to an infinite counter term. 



To cancel the pole in (6.9), we choose 



Hence, 



.(1) 



y 



(6.12) 



1pEFT 



d=A 



Fig.6(a) + Fig.6(b) 



1 + 



6 (0) 



^ + 7 + log7r + 21og| + 0(e) ) + ... 



. (6.13) 



On the other hand, expanding the full solution (4.16) in the vicinity of (i = 4, we 
obtain 



d=4 



1 + 



41og2-l p 
h 2 log — I + . . . 

4 po 



(6.14) 



First thing to note is that EFT correctly reproduces the leading log p term. Finally, match- 
ing the results, yields 



(L) = 



+ ^ log — 

L=po O Po 



AO) 



L = pQ 



^ A 2 7 



(6.15) 



The corresponding classical RG flow is thus given by 



/3. 



'A2 



L 



dX 



(0) 



dl' ~Y 



(6.16) 



All the above results are inherent in our master formula (4.17) as is evident from (5.6). 



Subleading corrections. When d = 2,3 there is no need to compute new diagrams. 
Indeed, in this case A2 ~ with 2 < a < 3, therefore only diagrams in figures 6(a) and 
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6(b) are relevant for our needs. Using (6.5) and (6.8), we obtain 



i^EFT 
i>EFT 



d=2 



d=3 



i+(^V+^^+... 

p J vr^ p6 



p J 8^2 p7 

Comparing with the asymptotic expansion of the fuh solution (4.16), yields 



As 
As 







d=2 



d=3 



87r2V3r2(5/3) 
15r(7/6)r(13/6) 



Po 



(6.17) 
(6.18) 

(6.19) 
(6.20) 



These results agree with (4.17). 

The last case to consider is d = 1. In this case, the full solution degenerates into a 
polynomial of fourth order in Z 



d=l 



(6.21) 



In particular, there is no need in any computation to fix As. Indeed, since as argued 
above, in isotropic coordinates the small parameter is m^, whereas according to (2.13) 
As ~ in d = 1, we conclude that there are no diagrams with mass insertions which scale 
in the same way as figure 6(a). Furthermore, the full solution does not contain odd powers 
of m ~ Po- As a result, we conclude that A3 must be zero to match the full solution in this 
case. Of course, this result fits our general formula (4.17). 

It is instructive to show how the EFT reproduces the full solution (6.21). This requires 
computation of one additional diagram apart from what we have computed so far, i.e. figure 
6(c) 



Fig.6(c) 



3{d - 2) 



4d 



.22) 



(d-4)(3d-4) 

For d = 1 this diagram reproduces the last term in (6.21). Note also that as expected it 
vanishes in d = 2. Indeed, in this case the full solution (6.21) truncates at second order in 
Z, therefore the overall contribution at order must vanish. As a last test of this result 
we checked that it produces correct contribution in d = 3. 



7 Summary of results 

Results. Our main result is the determination of the induced gravito-static polarization 
constants, or gravito-electric Love numbers, for any space-time dimension d, given by (4.23) 
together with the text below it and the normalization factor (3.5). 

We also computed the Love numbers for a free scalar field given essentially by (4.17). 
In d > 4 these describe also the tensor mode perturbations of the spatial metric. 
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Tests. The results passed two tests 

• The 4d results are consistent with [8, 9]. Their vanishing is dictated by the tanvr/ 
factor. 

• The residue of the pole at ^ = 1/2 was compared with the EFT thereby confirming 
the overall normalization factor. 

Main implications. 

• We settled a disagreement in the literature whether the calculation described in [8] 
and [9] implies a vanishing Love number or not. We find that in hindsight it indeed 
implies so and we supply the missing argument. This means that even though the 
horizons of 4d black holes do deform they are nevertheless "infinitely rigid" in this 
sense. 

• We answered negatively (at least for (i > 4) a conjecture about the absence of non- 
minimal world-line couplings in the BH effective action. 

• Our results for higher dimensions d > 4 are unquestionably novel, as well as those 
for the free scalar field. 

• We observed that the Love numbers are negative for certain ranges of parameters 
and interpreted that as indication for a novel non-spherical instability. 

See more in the discussion subsection 4.2. 



On the way we touched upon the following points of interest 

• We defined the static sector of the black hole effective action (2.2). 

• We defined / (3.7) and noticed that most of the analysis depends on /, d only through 
it. 

• We found quite a simple action (4.21) which encodes the full Zerilli potential. 

• We noted that there is an interesting limit to the theory as / — ?• oo where only the 
tanvr/ term survives in (4.23). 

• We provided a nice example of a classical RG flow, analyzed in sections 5,6. 
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A Matching normalization constant 

In this appendix we evaluate the normalization constant N (3.4). This is done by solving 
the EFT equation of motion in the presence of the A term (2.9) (obviously in flat space- 
time) and reading off the A coefficient defined in (3.2). 
The relevant part in the EFT action reads 




(A.l) 



where IS^'^/;! is defined in (2.9). The following analysis does not change by a multiplication 
of the action by an overall constant. The equation of motion is 




(A.2) 



Following (3.2) we write the solution as 




(A.3) 



Substituting back into (A.2) we find that 



N = {-y (coci) 



-1 



(A.4) 



where 



• Co is defined to be such that 




(A.5) 



namely, Gq 



co/r^ is a Green function satisfying AGq = 6{x). We have 



Co = 



1 



'd+l 




(A.6) 



where 0? 



'd+l ■ 



Vol ( S'^^ 
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ci is defined such that 

1 

rpd j'2l-\-d 

where I = (ii, ■ ■ ■ ,ii) and all the indices ii, . . . , are distinct. We find 

.. = (-2)'(|)^ = (-2y^ (Ai 



Substituting back into (A. 4) we obtain the requested quantity 



N = —/^ = ^ (A.9) 



B Useful formulae 

We collect some useful facts for the microscopic calculation, starting with the hypergeo- 
metric function, see for example [21]. 

The hypergeometric equation for u = u{y) is 

y{l -y)u' +[c- {a + b+ l)y]u' -abu = (B.l) 

The characteristic exponents are encoded in the Riemann P-symbol 

/ 



uiy) = P 



The hypergeometric series is defined by 




(B.2) 



(c)kkl 

{a)k :=^%^ = a-{a + l)...{a + k-l) (B.3) 
r(a) 



Pfaff's identity is 



F{a, b, c; X) = (1 - X^F (^c - a, b, c; (B.4) 

Expansion around X = is given by [21] (8.3.6) 

F{a,b,c;l-X) = ^^2lI^^-^^F{a,b, a + b + 1 - c; X) + 
1 (c — a)L [c — b) 

+ Mn^±i^X^--^ F{c-a,c-b,l+c-a-b;X) (B.5) 

r(a)r(6) 
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The Legendre polynomial is related to the hypergeometric function through 



Pl{x)=F{-l,l + 1,1; 



1 — X 



(B.6) 



2 



Finally we record some useful Gamma function identities 




(B.7) 
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